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Abstract 

We discuss the eigenvalue spacing statistics of the Glauber matrix for various 
models of statistical mechanics (a one dimensional Ising model, a two dimensional 
Ising model, a one dimensional model with a disordered ground state, and a SK 
model with and without a ferromagnetic bias). The dynamics of the one dimen- 
sional Ising model are integrable, and the eigenvalue spacing statistics are non- 
universal. In the other cases, the eigenvalue statistics in the high temperature 
regime are intermediate between Poisson and G.O.E (with P(0) of the order of 
0.5). In the intermediate temperature regime, the statistics are G.O.E.. In the low 
temperature regime, the statistics have a peak at s = 0. In the low temperature 
regime, and for disordered systems, the eigenvalues condense around integers, due 
to the fact that the local field on any spin never vanishes. This property is still 
valid for the Ising model on the Cayley tree, even if it is not disordered. We also 
study the spacing between the two largest eigenvalues as a function of temperature. 
This quantity seems to be sensitive to the existence of a broken symmetry phase. 
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Abstract 

Nous dicutons les statistiques d'ecart entre valeurs propres de la matrice de 
Glauber pour differents modeles de mecanique statistique (modele d'Ising unidi- 
mensionnel, modele d'Ising bidimensionnel, modele unidimensionnel avec un etat 
fondamental desordonne, modele SK avec ou sans biais ferromagnetique) . La dy- 
namique du modele d'Ising unidimensionnel est integrable, et la statistique d'ecart 
entre valeurs propres est non universelle. Dans les autres cas, la statistique de 
valeurs propres dans le regime haute temperature est intermediaire entre la loi de 
Poisson et la loi G.O.E. (avec P(0) de I'ordre de 0.5). Dans le regime intermediaire, 
les statistiques sont G.O.E.. Dans le regime de basses temperatures, les statis- 
tiques presentent un pic pour s = 0. A basse temperature et pour des systemes 
desordonnes, les valeurs propres condcnscnt autour des entiers, a cause du fait que 
le champ local sur aucun spin ne s'annule jamais. Cette propriete est encore vraie 
pour le modele d'Ising sur I'arbre de Cayley bien qu'il ne soit pas desordonne. Nous 
etudions egalement I'ecart entre les deux plus grandes valeurs propres en fonction 
de la temperature. Cette quantite semble etre sensible a I'existence d'une phase 
avec brisure de symetrie. 



1 Introduction 



The ideas of level spacing statistics emerged for the first time in the context of nuclear physics 
[0 0) where Wigner proposed computing statistical quantities from consideration of deter- 
ministic spectra. Later, these ideas were applied to quantum systems whose classical analogs 
are chaotic 0] [||. The idea of level spacing statistics is to calculate the difference between two 
consecutive levels, and to study the probability of occurrence P{s) of a level spacing s. The 
different generic behaviors of P{s) are classified according to random matrix theory [^] [^]. A 
generic case is the integrable spectrum. Each level is labeled by a set of quantum numbers, 
the energy levels are decorrelated and the statistics are Poissonian: P{s) = exp (— s). If the 
number of conserved quantities is too small, it is not possible to find a set of quantum numbers 
for each level, and the levels are correlated, that is, there exists level repulsion. The repulsion 
is linear and the level spacing statistics have the Gaussian Orthogonal Ensemble (G.O.E.) 
shape: 

Pis) = |.e-^^^ (1) 

If time reversal invariance is broken, and if the system is chaotic, the level spacing statistics 
has a Gaussian Unitary Ensemble (G.U.E.) shape: 

PM = 5^-^'. (2) 

where the repulsion is quadratic. The ideas of quantum chaos have been applied to various 
fields of condensed matter physics, such as disordered systems . Another field of application 
is strongly correlated electron systems [|] where the hope is to extract more information 
from finite size systems. In the present paper we wish to analyze the dynamics of classical 
spin systems using eigenvalue spacing statistics. We use here the term 'eigenvalue' rather than 
'level' since there are no energy levels as in quantum mechanics. We consider the 2^ x 2^ 
Glauber matrix, with the number of Ising spins, and diagonalize it for small clusters. The 
only symmetries are the lattice symmetries, that we treat using group theory, and the global 
Z2 symmetry. We can thus only study the dynamics for a small number of sites, typically on 
the order of 10 sites. The eigenstates are not physical, except for the Boltzmann distribution 
which corresponds to the upper eigenvalue A = 0. The other eigenvectors are not probability 
distributions, since the sum of their components is zero, so that their interpretation is not 
obvious. 

We have studied two different quantities. The first quantity is the distance A]\r{(3) between 
the two largest eigenvalues (one of them being zero for all (3) for A^ sites. In the infinite 
temperature limit, AAr(O) = 1 (see below for a proof of this fact). Notice that in what follows 
one unit of time corresponds to a single spin flip step, whereas a Monte Carlo Step would 
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correspond to N single spin flip steps. As a consequence, with the normalized time units 
the spectrum lies in the interval [—1, 1], and the normalized distance between the two largest 
eigenvalues is A]\f{(3) = An{(3)/N. In particular, AAr(O) = in the infinite temperature 
limit. One would expect the function An{P) to decrease as the inverse temperature increases 
since the relaxation times are expected to be larger for smaller temperatures. In the presence 
of a broken symmetry (for instance in the case of the two dimensional Ising model) , one expects 
that, in the thermodynamic limit, Aoo(/3) = if /3 > /3c since the broken symmetry states 
becomes degenerate with the Boltzmann distribution in the range of temperatures P > (3c- 
By contrast, in the absence of a broken symmetry (for instance for a one dimensional Ising 
model or for a model with a disordered ground state), the quantity Aoo(/3) should be finite, 
except in the limit (3 — +oo. However, our finite size study is far from the thermodynamic 
limit since we could only diagonalize the Glauber matrix for about ten sites. We were not 
able to distinguish between the different conjectured behaviors in the thermodynamic limit. 
Nonetheless, we study finite size effects and show that AAr(/5) decreases as the number of 
sites N increases. For some models (one dimensional Ising model and the frustrated one 
dimensional Ising model), we find that A]^{j3) is close to an exponential. In the case of the 
two dimensional Ising model, Ajv(/5) is clearly not exponential. 

Another quantity of interest is the eigenvalue spacing statistics of the full spectrum. We 
first consider the one dimensional Ising model. In this case, the dynamics are shown to be 
integrable. The corresponding eigenvalue spacing statistics are found to be non-universal, 
with a peak at s = which increases as the temperature decreases. In the case of the two 
dimensional Ising model with nearest neighbor coupling, the statistics are intermediate be- 
tween Poisson and G.O.E. for very small /3, with P(0) ~ 0.5. As /3 increases, the statistics 
evolves towards a G.O.E. shape (/3 ~ 1). At low temperatures no eigenvalue repulsion exists, 
and the statistics exhibit a peak at s = 0. The weight of the peak increases with /3. Next, we 
consider a frustrated one dimensional model with an extensive entropy at low temperatures. 
The evolution of the eigenvalue spacing statistics is similar to the case of the two dimensional 
Ising model. In the case of the Sherrington-Kirkpatrick (SK) model, we also have the fol- 
lowing evolution of the eigenvalues spacing statistics: no repulsion at very low temperatures, 
repulsion for intermediate temperature {(3 ~ 1), and no repulsion at low temperatures. An 
important property of disordered models is that their eigenvalues condense around integers 
at low temperatures. This is due to the fact that, except for disorder realizations with zero 
probability measure, the local field on any site is never zero. However, such a behavior also 
exists for non random systems, for instance for the nearest neighbor Ising model on the Cayley 
tree. 
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2 The Glauber matrix 



Glauber dynamics [|n[] are a single spin flip dynamics with a continuous time. If p({cr},t) is 
the probability to flnd the spin sytem in the conflguration {a} at time t, the master equation 
for the single spin flip dynamics is 



d / N \ N 

(3) 

The single spin flip transition probabilities are deflned as the probabilities that the spin aj 
flips from cTj to — (Tj while the other spins remain flxed. Since the Boltzmann distribution is a 
flxed point of the dynamics (|^), the transition probabilities have the form 

w,{{a}) = hl-a,tanh{f3J ^j))^ (4) 

where V{i) is the set of neighbors of the site i. Denoting the 2^ vector of the p({cr}, t) as p(t), 
equation (H) can be written as 

|p(t) = Gp(t), (5) 

where the matrix G is the Glauber matrix. Since the Boltzmann distribution is a steady state 
of the dynamics, its corresponding eigenvalue is zero regarless of temperature. The matrix G 
is not symmetric. It can however be related to a symmetric matrix M. To do so, we notice 
that the Glauber matrix satisfles the detailed balance, that is G^./jp^^"* = Gjj^aPa^ where p*^") 
is the Boltzmann distribution. As a consequence, 

(pi"')"" {pfY' = (pS?T"''g* (pL"')"'- (6) 
We call M the matrix deflned by 

M„,= {p<^>y"''G,,(pfy'\ (7) 
Then, M is symmetric. If p is a right eigenvector of the Glauber matrix, then 

J2 Ga/3P/3 = Ap^ (8) 

is equivalent to 

EM.,(p(°))-^/^, = A(pi°))-^/^., (9) 

SO that (p^^-*) ^ Pa is an eigenvector of M. We conclude that G is diagonalizable, and that 
all its eigenvalues are real. 
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The spectrum in the infinite temperature hmit can be understood as follows. If we call 



IV') = E/(W)ki) 



fat), 



then the dynamics reads 



N 



i=l 



dt 



so that the eigenvalues of the Glauber matrix at infinite temperature are of the form 



N 



N 1 ^ 

1=1 



(10) 



(11) 



(12) 



where //^ = ±1. The spectrum in the infinite temperature limit is thus made of eigenvalues at 
integer values between — and 0, with a degeneracy given by the binomial coefficients. 

Another property of G is that for bipartite lattices, such as the square lattice or the 
Cayley tree, the spectrum of G is symmetric: if A belongs to the spectrum then — — A is an 
eigenvalue also. The proof is as follows. Let Xjcr} be an eigenvector of M, with an eigenvalue 
A: 



N 



1 ^ 1 

AX{(7} = - ^ -(1 - a, tanh (/3 J/i,))X{(7} + ^ ^ T^'^^'^'' ^^^^ 

1=1 i=l r-^ * 



where hi is defined by 



Let Y{cr} be defined as 



Y{<7} = (-l)^Mx{a}, 



(14) 



(15) 



where v{o^ is the number of up spins in the configuration {cr}. {cr} is deduced from {cr} by 
fiipping the spins of one of the two sublattices. Then, 



N 1 



(16) 



1=1 

N 



_i_ y^(^_2^t/{o-i,...,-o-i, 



,crN} _ 



i=l 



2 cosh (pJhi) 



X{a-i, —ai, ctm} 



= (-l)-W 



TV 



N 



-E 



J2-{l + aite.nh{/3Jhi))-K{a} 

i=i ^ 



-X{a-i, —Gi, ...jO-jv} 



(17) 



^ 2 cosh {(3Jhi 
{N + A)(-l)'^^'^>X{a} = -(iV + A)Y{a}. 



(18) 
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Given an eigenvector X for the eigenvalue A, we have built an eigenvector Y for the eigenvalue 
-A^- A. 

The difference between (^) and the Schrodinger equation is that quantum mechanics pre- 
serves the scalar product, leading to Hermitian Hamiltonians. More over, the space of physical 
states is a Hilbert space, and each state of the Hilbert state is physical. In the case of the 
Glauber matrix, no vector space is present in the sense that the sum of two probability dis- 
tributions is not a probability distribution. However, some quantities are conserved by the 
dynamics. It is easy to show that the eigenvectors of G for the non-zero eigenvalues have the 
property that 

EpW = 0- (19) 
This is a simple consequence of the fact that the Glauber matrix preserves the quantity 

EpW- (20) 

3 One dimensional Ising model 
3.1 Integrability of the dynamics 

In the case of the one dimensional model, the Glauber dynamics is integrable. To show this, we 
follow Glauber and write evolution equations for the correlation functions. We call {t) 
the n-point correlation function, with ia ^ ip ii a ^ (3, that is 

Rff...,^.it) = {E^^At)). (21) 

a 

Then following Glauber, we write the evolution equation of -Rj-"^ i^{t) under the form 

Jt^nL^it) = -2(a,(t)...a,„(t) ^({a}) + ... + ^,„({a}))), (22) 

where the transition probabilities are given by (^). We notice that in the one dimensional 
case, each spin has two neighbors, so that Wi{{a}) can be written as 

w.iW}) = ^ (l - ^^.(^.+1 + ^.-i)) , (23) 

with 7 = tanh 2/3 J. Here, we take periodic boundary conditions, but the case of open boundary 
conditions is similar. Inserting ( P5D into (^) we get 

J n 

j,i^tAt) = -nitUi^ + i 5: E (^^.+. n ^^.)- (24) 
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The terms with e = 1 collect the right neighbors, and e = — 1 corresponds to the left neighbors. 



The correlation function in leads to a (n-2)-point correlator if 3(3, i^+t = "ifs or to a n-point 



correlator if not. The expression (^4) can be brought under the form 



^ ( ) ( ) n I n ^ ^ 




If none of the sites ii,...,in are neighbors, the term containing vanishes. However, 

if at least two sites in the set ii, ...,in are neighbours, we have to take into account a term 
containing in the evolution of R^"'\ It is clear that ( ^5]) is nothing but a rewritting of 

(1^) in the case where all the sites have only two neighbors. The number of distinct correlation 
functions is 

^:^^l=2^ (26) 



n=0 



n 



which is equal to the number of spin configurations. The system ( P3| ) is integrable. Glauber 
gives the explicit solution for R[^\t). The equation giving dR^'^\t)/dt contains only linear 
combinations of R^'^\ In order to solve for the three points correlation functions, we inject 
Glauber's solution into the evolution equation for R^^\ diagonalize the associated matrix and 
get a first order differential equation, which is explicitly integrable and yields the second order 
correlation functions. The entire hierarchy can be solved by this method since dR^^^dt does 
not contain R^'^^ with p > k. 



3.2 Eigenvalue spacing statistics 

In order to calculate the eigenvalue spacing statistics, we need to take all the symmetries of 
the lattice into account. Here, the symmetries are so obvious that we do not require a group 
theory treatment. We work with an open Ising chain. This graph is invariant under the 
reflection and the identity operators. We denote the basis of our 'Hilbert' space as {cr}). We 
use quotes since there is no vector space structure on the probability distributions. However, 
to diagonalize the Glauber matrix, we can use an analogy to quantum mechanics. If R is the 
refiection operator, we form the combinations 

IW)e = ^(IW) + ei?|M)). (27) 

This operation leads to states with a well defined behavior under the refiection. The resulting 
state is either symmetric (e = 1) or antisymmetric (e = — 1). The antisymmetric state may be 
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zero if is invariant under the reflection. The dimension of the antisymmetric sector is 



where [] denotes the integer part. Finally, we take into account the global Z2 symmetry 
of the Glauber matrix. [The matrix elements of the Glauber matrix are invariant under 
the transformation {cTj} {— cTj}]. Taking into account all the symmetries, we diagonalize 
the Glauber matrix in the reflection symmetric and the antisymmetric sectors, and in the 
Z2 symmetric and antusymmetric sectors. The evolution of the eigenvalues of one sector as a 
function of the inverse temperature is plotted in figure |T] (a). No avoided crossings are present, 
which is what is expected for an integrable system. The difference between the two largest 
eigenvalues is plotted on figure |^ as a function of the inverse temperature. Eventhough some 
deviations are visible, the behavior of A7v(/3) is close to an exponential decay. Finite size 
effects are visible: A7v(/5) decreases as the number of sites increases. 

The eigenvalue spacing statistics are found to be non universal. For instance, -P(s) is 
plotted on figure | for /3 = l.The height of the peak at s = decreases as the temperature 
decreases. Eventhough the dynamics is integrable, the eigenvalue spacing statistics are not 
of the Poisson form (|l]). Such a behavior for the spectrum of integrable systems has already 
been found in the context of integrable quantum fluids (see 0). 

4 Bidimensional Ising model 
4.1 Dynamics 

In this case, the dynamics are no longer integrable. The evolution of the correlation function 
is not a linear equation, as it was in the case of the Ising chain. This is essentially due to the 
fact that, with four neighbors, one has to introduce a cubic term in iyj({cr}) given by equation 




(28) 



and the dimension of the symmetric sector is 




(29) 



tanh {(3J{ai + (J2 + 0-3 + 0-4)) = a{ai + (J2 + 0-3 + 0-4) + a'{ai + 0-2 + 0-3 + (J4) 



3 



(30) 



with 



a' 



a 



— (8 tanh 2/3 J - tanh 4/5 J) 
(tanh 4/3 J - 2 tanh 2/5 J) . 

48 



(31) 



(32) 
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In the one dimensional case, we could integrate the dynamics because dR^^ydt was only a 
function of i?^^-* with k < n. In the two dimensional case dR^'^^dt is also a function of R^'^^ 
with k > n, so that the hierarchy is no longer integrable by this method. It is not because 
one does not know how to solve the dynamics that these dynamics are not integrable. The 
analysis of the spectral statistics of the Glauber matrix may be useful to determine whether 
or not there exists some conserved quantities in the Glauber dynamics of the the Ising model. 



4.2 Use of group theory 

We use group theory to find the symmetries of the clusters for which we shall diagonalize the 
Glauber matrix. Notice that we are restricted to small sizes since the size of the 'Hilbert' 
space is equal to 2^. In practice, and to have reasonable execution times, we are restricted 
to < 13. The first step is to determine the symmetry group of the lattice, that is to 
enumerate all the permutations that leave the lattice invariant. To do so, we do not test all 
the A^! possible permutations since the computation time may be huge. Instead, we use the 
following procedure. We first label the lattice sites and give the list of bonds. We determine 
all the possible images a{l) of the site 1, that is all the — 1 sites. Then, for each of the 
possible images of the site 1, we determine the images a{2) of site 2 which respect the lattice 
symmetry: if there is a bond between 1 and 2, there must be a bond between cr(l) and cr(2). 
If there is no bond between 1 and 2, there must be no bond between a{l) and cr{2). At 
this point, we have a list of potential permutations beginning with a{l) and o"(2). Then, 
we determine all the possible images of site 3 which leave the lattice invariant. We thus get 
a tree structure, but, during the construction, some branches shall stop. At the end of the 
process, that is when a{N) has been determined, we get all the permutations which leave the 
graph invariant. The second step is to determine the classes and the table of characters of the 
group. We use the program in [Q, which automatically determines the classes and the table 



of characters. In a third step, we have to determine the size of the blocs corresponding to the 
irreducible representations, and how many times a given irreducible representation appears. 
The dimension of the blocs corresponding to the irreducible representation (j) is equal to 

Dtm^^^ = Wx{9)x'^'\~9). (33) 
^ sec 

g is the representation of the group element g in the 'Hilbert' space, h is the cardinal of 
the group G, xid) is the trace of g and x^^Kd) is read from the table of characters at the 
intersection of the line corresponding to the representation (j) and the column of the class of 
g. The fourth step is to implement a Gram-Schmidt procedure to determine the basis of one 
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block corresponding to the representation (j). We first use the projector 

pO) = ^^0)(^)^. (34) 

geG 

A basis element of the Hilbert space is coded as a binary number of size A^. Zero corresponds 
to a down spin, and 1 codes an up spin. In order to label the basis vectors, we use the decimal 
representation of the binary number of size A^. We denote the corresponding vector \ipk)- The 
procedure consists of scanning all the states \tpk) and to determine ko such as P^^'^\4)k) = 
ii k < ko and P^^''\ipko) 0. The state P^^^lipko) is the first vector of the basis that we are 
looking for. Once we have found the first vector of the basis, we continue to scan all the states 
\ipk), but we project them with 

Po^ = T.{^ko\mkoy9- (35) 

g€G 

If PQ^lipk) = 0, we forget about \ipk) and project {ipk+i)- If Po^l'^k) 7^ 0, we try to incorporate 
Po^li^k) into the basis using a modified Gram-Schmidt procedure |jl3|]. If PQ^'l'ipk) is a linear 
combination of the basis vectors, then we discard it and project iV'fc+i)- If it is not, we 
incorporate it into the basis, after having orthogonalized it, and we make the projection test 



for |'?/'fc+i). At the end of the procedure, the dimension of the basis must be equal to (|33D. We 
note that it is not possible to store all the components of the orthonormal basis because of 
limited storage capacity. In order to save memory, we stored only the non zero components. 
The fifth step is to take into account the global Z2 symmetry. The sixth and last step is to 
diagonalize the Glauber matrix using the basis that has been determined at the fourth step. 
The size of the matrices to be diagonalized are small enough, so that we can use the Jacobi 
method. 

4.3 Results 

We work with a 3x4 lattice with periodic boundary conditions. The number of representations 
is equal to 15, and the maximal block dimension is 335. The spectrum in a given sector of 
symmetry of the 3x4 square lattice is pictured in Figure |I| (b) as a function of the inverse 
temperature. In the limit (3 —>■ 0, we recover degeneracies for integer eigenvalues (see relation 
(O)). As the inverse temperature increases, the degeneracies present for /3 = are lifted. 



but the eigenvalues from two different degeneracies are not free to cross, due to eigenvalue 
repulsion. 

We studied the evolution of A^iP) as a function of /5 for the 3x4 lattice and for the 3x3 
lattice. The result is plotted in figure Finite size effects are visible; A^r(/?) decreases if 
A^ increases, and the evolution of A]^{/3) is clearly non exponential. In the thermodynamic 
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limit, Aoo(/5) = if P > Pc and Aoo{P) > if P < Pc, which means that lnAoo(/5) — > -oo if 
P ^ P~ . Our finite size study is consistent with such a behavior. However, it would be useful 
to analyze larger samples, which will be done in the near future. 

We now discuss the shape of the eigenvalue spacing statistics P{s). If the inverse tem- 
perature is very small, the degeneracies of the P = case are lifted and the once degener- 
ate eigenvalues spread out linearly. The corresponding statistics are plotted in Figure |^ for 
P = 0.01. In this case, we find -P(O) ~ 0.5 and the statistics are close to the Poisson law for 
large s. If P increases, one reaches the eigenvalue repulsion regime (see figure ^ (b)). Linear 
eigenvalue repulsion is visible on figure ^ since, for /? = 1, P(s) ~ s for small s. However, for 
s of order unity, large deviations to the G.O.E. law occur. For large P, the statistics are not 
universal (see Figure |^ for P = 10), with a peak at s = 0. The weight of this peak increases 
as P increases. 



5 Frustrated one dimensional model 
5.1 The model 

We consider the one dimensional antiferromagnetic Ising model with antiferromagnetic next- 
nearest-neighbor interactions. This model can be seen as a succession of triangles, as pictured 
in Figure ^ and can be solved via a transfer matrix formalism, with the sites gathered as shown 
on figure |^. The transfer matrix has the form 



t=[b a)' P6) 

with 

A - 1 j ,B-y ^ j . (37) 

In (p6|), the states are ordered in the form | t)T)) I T)i)) I I lA)- Because of the 

form (^) of the transfer matrix, if (?/), if) is and eigenvector of T, then {ip + ip^ip + if) and 
{ip — <f,ip — <f) are eigenvectors of T for the same eigenvalue, so that the eigenvalues of T 
are the eigenvalues of A + B and A — B, and the initial 4x4 problem is reduced to two 2x2 
problems, due to the time reversal invariance. The partition function is simply Zjy = TrT^^"^ 
for a site chain (A^ is assumed to be even). The zero temperature entropy is found to be 
extensive, of the form S{Q)/N = In 2/2, whereas in the corresponding ferromagnetic problem, 
the entropy is finite at zero temperatures. This one dimensional antiferromagnetic model has 
thus the same properties as the triangular antiferromagnet |]14|, namely the number of zero 



temperature ground states is proportional to expaA^, with a a constant. 
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5.2 Results 



We work with an open chain version, so that the only symmetries are the inversion and the 
global Z2 symmetry. We have already explained how to treat these symmetries (section |3.2| ). 

The evolution of A7v(/5) as a function of the inverse temperature /3 is plotted in Figure 
|] for different values of A^. We observe finite size effects: Ajv(/5) decreases with A^. The 
variations of AAr(/?) as a function of /5 are linear in a semi-log plot, so that A7v(/5) decreases 
exponentially with /5. 

The evolution of the eigenvalues of the Glauber matrix as a function of the inverse tem- 
perature is plotted in figure |I| (c). Level repulsion is visible for 0.2 < /3 < 1. 

The eigenvalue spacing statistics in the high temperature regime is plotted in Figure ^ 
(/3 = 0.01). In this range of temperatures, the statistics are intermediate between Poisson 
and G.O.E. P(s) has a maximum, but is occurs for smaller separations than in the G.O.E. 
case. If the inverse temperature increases, the maximum occurs for larger separations, and 
the eigenvalue spacing statistics is close to the G.O.E. law for /3 = 1 (see Figure]^). For larger 
temperature (/? = 3 on Figure |l^), -P(O) is not zero (for instance, -P(O) ~ 0.4 for /? = 3) and 
P(s) is close to the Poisson law for large spacings. 



6 SK model 
6.1 The model 

This model was proposed in 1975 as an 'exactly solvable' spin glass model |]T5[. For general 



reviews on the problems of spin glasses, we refer the reader to references [T^ fl^ [|l^. The 



SK model is defined by the disordered, infinite range interaction Hamiltonian 

if = ^ JijOiOj, (38) 

with quenched random interactions Jij with a gaussian distribution 

In this section, we study the symmetric case 5 = 0. The case 5 > 1 is the object of section ^ 
We refer the reader to the reviews previously quoted for the solution of this model. For T < J, 
the model is a spin glass, and a paramagnet for T > J . In the glass phase, the model has 
a large number of thermodynamic phases, with no symmetry connecting ground states. The 
different ground states are separated by large barriers, proportional to the system size. An 
other feature of glassiness is the presence of aging, associated with slow relaxation processess. 
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6.2 Eigenvalue spacing statistics 

In the case of the SK model, we do not need to look for the lattice symmetries since, except 
for some special cases of zero probability, the random infinite range interactions break all the 
lattice symmetries. The only symmetry to be taken into account is the global Z2 symmetry. 

The evolution of the eigenvalues as a function of the inverse temperature is plotted in 
Figure |Tl| (a) in a given symmetry sector. We observe that, in the limit (5 — > +00, the 
eigenvalues condense around integers. This property is due to the fact that, in this limit, the 
symmetric representation M, of G, is diagonal since the local field never vanishes (except for 
some disorder configurations of zero probability measure), and that the diagonal coefficients 
are integers. This cluster property in the zero temperature limit is an effect of disorder, and 
is not related to the existence of glassiness. However, as we shall see in section ^, we find a 
similar behavior for a glassy ferromagnetic system with no disorder ! 



The evolution of A7v(/3) as a function of /3 is plotted in Figure |T2[ The variations of 



AAr(/5) as a function of /3 are approximately linear in a semi-log plot, so that AAr(/3) decays 
exponentially as a function of /?. 

The eigenvalue spacing statistics in the high temperature regime are intermediate 

between the Poisson law and the G.O.E. law (see Figure ^ for /5 = 0.01). In this regime, 
P(0) ~ 0.5 and -P(s) is close to the Poisson law for large separations. If the inverse temperature 



/? increases, the maximum of -P(s) occurs for larger separations (see Figure |T3| for /5 = 0.1) 



and finally -P(s) reaches the G.O.E. law (see figure for /5 = 1). In the very low temperature 
limit, the statistics are close to the Poisson law, with a peak for s = (see Figure |1^ for 
/3 = 10). The height of this peak increases as /? increases. 

6.3 Zero temperature density of zero eigenvalues 

At zero temperature, the eigenvalues of the Glauber matrix of the SK model are integers, and 
the Glauber matrix is diagonal in the natural basis, excepted for some disorder configurations 
of zero probability measure. As a matter of fact, some eigenvalues are zero. The states 
corresponding to the zero temperature zero eigenvalue are the metastable states of the SK 
model, and their degeneracy is expected to increase exponentially with the system size, as 
expected from the number of solutions to the TAP equations [0. It is thus interesting 
to calculate the density of zero eigenvalues in the zero temperature limit as a function of the 
system size, even if our approach is restricted only to small system sizes. First, we analytically 
treat the case of the 3 site SK model, as a warm up exercise. 
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6.3.1 Warmup: the 3 site SK model 

We call (Tj the spin variable on site i (i = 0, 1, 2), and J {K,L) the random bond between ctq 
and o"! ((Jo and a2, <7i and (T2). Since the local field never vanishes, the Glauber matrix is 
diagonal, and the diagonal coefficients are 

Ai = A{+,+,+) = A{-,--) = -^ + ^e{J + K) + ^e{J + L) + ^e{K + L) (40) 

A2 = A(+,+,-) = A(-,-,+) = -^ + ie(J-X) + le(J-L)-ie(ir + L) (41) 

A3 = A(+,-,+) = A(-,+,-) = -^ + ie(i^-J)-ie(J + L) + ie(i^-L) (42) 

A4 = A(-,+,+) = A(+,-,-) = -^-l6(J + X)-^6(J-L)-l6(K-L), (43) 
where e{x) denotes the sign of x. The gaussian bond distribution is replaced by 

Ps'iJ) = + S'- J)9{J +^-S') (44) 

in order to analytically perform the forthcoming integrations. The variable S' is equal to 
5/3. After lengthy but straightforward calculations, we get the probabihty Vi{6') that Aj 
(i = 0, 1, 2, 3) takes the value zero: 

• ii5'> 1/V3, Vo{5') = 1 and Vi{5') = 0. 

• if < 5' < 1/V3, 

Vo{5') = \+'-^S' + l6--^-^5- (45) 
V^i5') = (46) 

• if-l/V3<(5'<0, 

V.{5') = (48) 

• finally, if 5' < -1^3, Vo{5') = and Vi{5') = 1/3. 

Moreover, Vi = V2 = V3. We observe that Vo + 3Vi = 1 regardless of the value of the 
ferromagnetic bias 6'. In the ferromagnetic region, only one state (ferromagnetic) has a zero 
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eigenvalue. In the frustrated region, three 'disordered' configurations share the zero eigenvalue, 
while the probability of the ferromagnetic state is vanishing. For this three site toy model, 
the probability of finding a zero eigenvalue is a constant as a function of the ferromagnetic 
bias. This example shows that two mechanisms are responsible for the appearance of zero 
eigenvalues: (frustration + disorder) or broken symmetry. Notice that both frustration and 
disorder are important here since without disorder, the Glauber matrix would not be diagonal 
in the zero temperature limit. 

6.4 N site SK model 

We studied the total probability Pat (5) of finding a zero eigenvalue at zero temperature. This 
quantity is plotted in Figure |l^ in the frustrated region (5 = —2), in the zero bias case (5 = 0) 
and in the ferromagnetic region {6 = 2). We observe that in both cases Vn{6) increases with 
N, but more slowly in the ferromagnetic case than in the zero bias case. This suggests that 
the (frustration + disorder) mechanism to generate zero eigenvalues of the Glauber matrix 
at zero temperature is more efficient than the broken symmetry mechanism, in the limit of 
large N. Moreover, we see that the variations of Vn{6) as a function of are compatible 
with an exponential growth. However, we are restricted to small systems, so that it would be 
interesting to do the same study for larger systems. This shall be done in the near future. 
Another quantity of interest is the complete distribution function P{An{i3)). This quantity 
is plotted in Figure [1^ in the low temperature and high temperature regime. We observe 
that, in the low temperature regime, the standard deviation of Ai^{f3) is of the same order as 
AnI/S), whereas in the high temperature regime, the standard deviation is small compared to 
the mean value. 

7 SK model with a ferromagnetic bias 

We now consider the SK model with a ferromagnetic bias 6 > 1. In this regime, the SK model 
possesses a paramagnetic/ferromagnetic transition as a function of temperature. The spectrum 
of the Glauber matrix is plotted on Figure |ll] (b) as a function of the inverse temperature 
/3. Level repulsion is visible, and the eigenvalues are attracted by integers in the limit of 
large /3, which is consistent with the fact that the Hamiltonian is disordered. We plotted on 
Figure |12| Aj^{j3) for different values of the bias 6. We observe that A]\f{P) decreases as the 
ferromagnetic bias increases. This is due to the fact that the paramagnetic/ferromagnetic 
transition temperature of the biased SK model increases with the bias 6 |]T5[, so that the 
inverse transition temperature decreases as 6 increases, leading to a decrease of A]\f(P) as 5 
increases. Moreover, the variations of Ajv(/9) become non-exponential if the ferromagnetic bias 
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increases (see the case 5 = 3 of Figure O), which is consistent with what has been observed 
in the case of the two dimensional Ising model. The eigenvalue spacing statistics are plotted 
in Figure 17 and Figure IS. Level repulsion appears even for small values of (3. 



8 Ferromagnetic Cayley tree 

Trees were introduced in statistical physics as early as the 1930's in order to implement mean 
field theories [1^ |2^. In this case, only the central spin was considered and the border of 



the tree (which represents a finite fraction of the spins) was sent to infinity. This limit is 
the Bethe-Peierls limit. However, it is possible to develop statistical physics on a Cayley 
tree, with the border included. This was done in the seventies, and a continuous transition 
was discovered |21| |22|. The dynamics was studied only recently |]10|. A cross-over (not a 



transition like in the SK model) to a glassy regime was then found below the temperature 
scale Tg ~ J/ Inn, with n the number of generations. The temperature scale Tg decreases in 
the limit n +oo, but very slowly, so that the glassy regime exists even in the macroscopic 
regime. The zero temperature barriers scale like Jinn, with n the number of generations. 
If the coordination of the bulk sites is equal to 3, the local field of the bulk sites cannot 
vanish. This property is reminiscent on the behavior of disordered systems. The analysis of 
the spectrum of the Glauber matrix at low temperatures reveals that the clustering property 
also exists for the Cayley tree, as shown in Figure |1^, where we plotted the eigenvalues of 
the Glauber matrix on the Cayley tree at low temperature. This indicates a striking analogy 
between the ferromagnetic Ising model on the Cayley tree and disordered systems. 



9 Conclusion 

Since the contents of this paper were already summarized in the introduction, we end here 
with some final remarks and open questions. The one dimensional case is very special due to 
the existence of an underlying integrable dynamics. In the other cases, we have the following 
scenario for the evolution of the eigenvalue spacing statistics as a function of (3: there exists 
an intermediate regime (/? ~ 1) where the statistics are dominated by the non integrability of 
the dynamics, and are of the G.O.E. type. In the low temperature limit and in the high tem- 
perature limit, no eigenvalue repulsion was found. Another important feature of the spectrum 
of the Glauber matrix is that for disordered systems, the eigenvalues are close to integers in 
the zero temperature limit. This is a consequence of the fact that the local field never vanishes 
on any site. This property is not restricted to disordered systems, since we also observe eigen- 
value clustering at low temperatures for the nearest neighbor Ising model on the ferromagnetic 
Cayley tree. Notice that the clustering property is not related to glassiness since disordered 
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systems without glassiness may exist (the biased SK model with a sufficient ferromagnetic 
bias) and ferromagnets exist which are known to be glassy at low temperatures and which are 
not expected to have this property (regular fractals, percolation clusters, hyperboUc lattices). 

We also studied the separation between the two largest eigenvalues. This quantity exhibits 
strong finite size effects. In the case of the one dimensional Ising model, the frustrated one 
dimensional model, and the unbiased SK model, this quantity decays approximately exponen- 
tially with the inverse temperature. However, in the case of the two dimensional Ising model 
and in the case of a biased SK model, the decay is faster than exponential. This suggests that 
the separation between the two largest eigenvalues is sensitive to the existence of a second 
order transition. In the thermodynamic limit, Aiv(/3) is expected to diverge at the transition 
temperature in a semi- log plot. 

We also studied the density of zero eigenvalues in the SK model. At zero temperature, this 
quantity increases with the number of sites. We found variations that are compatible with an 
exponential growth. However, it would be interesting to investigate larger systems. As far as 
the highest non-zero eigenvalue probability distribution is considered, we found that at low 
temperature, the standard deviation of this distribution is of the order of the average value 
of the distribution, whereas in the high temperature regime, the standard deviation is small 
compared to the average value. This result is compatible with the fact that glassiness is due 
to an accumulation of eigenvalues in the vicinity of zero. 

In the present study, we only treated the case of small systems. A generalization to larger 
systems is under way, in which it would be especially interesting to study the separation 
between the two largest eigenvalues, as well as the probability of occurence of zero eigenvalues 
at zero (and non zero) temperatures. Another important point is that it would be interesting 
to have analytic expressions for Ajv(/3). This question is also under study. We shall also 
address in the near future the case of a Glauber dynamics with an external magnetic field. 

I thank B. Dougot for a critical reading of the manuscript, and I am indebted to J.C. Angles 
d'Auriac and H. Bruus for help with group theory. I also thank C. Denniston for reading the 
manuscript. I thank the referee for pertinent criticism. 
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Figure captions 



Figure 1: 

Evolution of the eigenvalues of the Glauber matrix as a function of the inverse temperature for: 
(a) an 8 site Ising chain in the reflection antiperiodic sector and the Z2 antiperiodic sector. No 

avoided crossings are present, as expected for an integrable system; (b) the 3x4 Ising model 
on a square lattice, in the representation number 3 and the antisymmetric sector. The 
evolution of 32 eigenvalues is plotted on this figure; (c) the frustrated one dimensional model 
with nearest neighbor interactions, for 8 sites in the antisymmetric reflection sector, and the 
antisymmetric Zi sector. 



Figure 2: 

Evolution of the parameter Ajv(/3) as a function of the inverse temperature /3, for different 
sizes (open Ising chain). 



Figure 3: 

Eigenvalue spacing statistics of the one dimensional Ising model with nearest neighbors cou- 
plings. The inverse temperature is /3 = 1. The statistics are not universal. 



Figure 4: 

Evolution of A^(/3) as a function of the inverse temperature /3 for the 3x4 lattice {N — 12) 
and the 3x3 lattice (A^ = 9). Finite size effects are visible, and the variations of Ajv(/3) are 
not exponential. 



Figure 5: 

Eigenvalue spacing statistics of the 3x4 square lattice for /3 = 0.01 and /3 = 1. The Poisson 
and G.O.E. distribution are plotted in dashed lines. In the regime where initially degenerate 

eigenvalues spread linearly (J5 = 0.01), -P(O) ~ 0.5 and the statistics are close to the Poisson 
law for large s. In the eigenvalue repulsion regime {j3 = 1), the statistics exhibit eigenvalue 
repulsion (P(0) = 0), but the shape of P{s) is distinctly different from the G.O.E. shape. 



Figure 6: 

Eigenvalue spacing statistics of the 3x4 square lattice for /3 = 10. The Poisson distribution is 
plotted in dashed lines. 
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Figure 7: 

The one dimensional Ising model with next-nearest-neighbor coupling, and its representation 
as a succession of triangles. The dashed lines represent how the sites are gathered in the 
transfer matrix formalism. 

Figure 8: 

Evolution of Ajv(/5) as a function of P for N = 6, 7, 8, 9, 10 in the case of the disordered one 
dimensional model of Figure |^ If > A^', < Ajy^P); An{P) decreases exponentially 

with the inverse temperature j3. 

Figure 9: 

Eigenvalue spacing statistics of the frustrated one dimensional model with next nearest neigh- 
bor interactions, for an inverse temperature f3 = 0.01 and (3 = 1. 

Figure 10: 

Eigenvalue spacing statistics of the frustrated one dimensional model with next nearest neigh- 
bor interactions, for an inverse temperature (3 = 3. 

Figure 11: 

Evolution of the eigenvalues as a function of the inverse temperature (3 (a) in the Z2 anti- 
symmetric sector of a 7 site SK model, and for a given disorder configuration; (b) in the Z2 
antisymmetric sector of the 7 site SK model with a ferromagnetic bias, and for a given disorder 
configuration. The ferromagnetic bias is S = 1.5. 

Figure 12: 

Evolution of Ajv(/3) as a function of (3 for = 9 in the case of the SK model, for different 
values of the ferromagnetic bias 6. An{(3) has been averaged over 200 realizations of the 
disorder. We plotted for S = 0,1, 1.5,2,3. An{P) decreases if the ferromagnetic bias 

S increases, which is consistent with the fact that the transition temperature increases as a 
function of 6. 

Figure 13: 

Eigenvalue spacing statistics of the SK model for an inverse temperature (3 = .01 and P = 0.1. 
The G.O.E. shape and the Poisson shape are plotted in dotted lines. 

Figure 14: 

Eigenvalue spacing statistics of the SK model for an inverse temperature (3 = 1 and (3 = 10. 
The G.O.E. shape and the Poisson shape are plotted in dotted lines. 

Figure 15: 

Probability to find a zero eigenvalue at zero temperature, as a function of the number of sites, 
for 5 = —2,5 = and 5 = 2. This is a semi-log plot. 
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Figure 16: 

Distribution of the largest non-zero eigenvalue — Ajv(/3) for a 9 site SK model, and for (3 = 0.2 
and (3 = 1.7. At low temperatures, the standard deviation is of the same order of magnitude 
as the mean value, whereas this is not the cas at high temperatures. 



Figure 17: 

Eigenvalue spacing statistics of the biased SK model for an inverse temperature = .01 and 
(3 = 0.1. The bias is 5 = 1.5. The G.O.E. shape and the Poisson shape are plotted in dotted 
hnes. 



Figure 18: 

Eigenvalue spacing statistics of the biased SK model for an inverse temperature (3 = 1. The 
bias is 6 = 1.5. The G.O.E. shape and the Poisson shape are plotted in dotted lines. 
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Figure 19: 

Eigenvalues of the Glauber matrix on the Cayley tree with 2 generations and coordination 3 
(10 sites). The inverse temperature is (3 = 2. 
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